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SATELLITE OVER A SELECTED EARTH POSITION 
By T. H. Skopinski and Katherine G. Johnson 


SUMMARY 


Expressions are presented for relating the satellite position in 
the orbital plane with the projected latitude and longitude on a rotating 
earth surface. An expression is also presented for determining the azi- 
muth angle at a given burnout position on the basis of a selected passage 
position on the earth's surface. 

Examples are presented of a satellite launched eastward and one 
launched westward, each passing over a selected position sometime after 
having completed three orbits. Incremental changes from the desired 
latitude and longitude due to the earth's oblateness are included in the 
iteration for obtaining the azimuth angles of the two examples. The 
results for both cases are then compared with those obtained from a com- 
puting program using an oblate rotating earth. Changes from the selected 
latitude and longitude resulting from incremental changes from the burn- 
out azimuth angle and latitude are also analyzed. 


INTRODUCTION 


In the recovery of an artificial earth satellite it is necessary to 
bring the satellite over a preselected point above the earth from which 
the reentry is to be initiated. It is the purpose of the present paper 
to determine the azimuth angle which must exist at a given burnout point 
In order to have the satellite vehicle pass over any particular earth 
position after a selected small number of orbital passes. For this pur- 
pose equations are developed which give the latitude and longitude of a 
point on the surface of a rotating spherical earth directly under the 
satellite at any time. An expression is developed which, with the orbit 
characteristics and the selected passage position on the earth’s surface, 
will establish after only a few Iterations the required azimuth angle at 
the burnout point. 

Elliptical orbits about a rotating spherical earth are first assumed. 
Although this simplification ignores the effects of the earth’s oblateness, 
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it is convenient to examine the oblateness effects later as incremental 
changes which can be accounted for in an iteration procedure. Departures * 

from the desired latitude and longitude passage position resulting from 
incremental changes in the azimuth angle and latitude at burnout are also 
analyzed. 

SYMBOLS 

semimajor axis of elliptic orbit 

eccentric anomaly, defined by equation (9) 

eccentricity of orbit, defined by equation (6) 

gravitational constant at earth surface 

inclination angle of orbital plane 

nodal point at which satellite crosses equator 

number of completed orbital passes, referenced to launch 
latitude 

0 center of earth 

P perigee location in orbital plane 

p semi 1 at us rectum of ellipse, defined by equation (l) 

R radius of earth 

r distance of satellite from earth center 

S position of satellite in orbital plane 

T satellite period, defined by equation (7) 

t(9) time from perigee, min 

V satellite velocity, defined by equation (3) 

V c satellite velocity for circular orbit, defined by equation (4) 

elevation angle between local horizon and velocity vector, 
positive upwards 
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0 angle In orbital plane measured from perigee pointj also, true 

anomaly 

A longitude measured east of prime meridian at Greenwich 

^NS = *N ~ ^S 
Mri = A n - Ap 

0 geocentric latitude measured north or south of equator 

i|r azimuth angle measured clockwise from North 

ft angle defining rotation of orbital plane about earth's axis 

(u angle in orbital plane between nodal point and perigee and, 

also, argument of perigee, deg 

ajg angular velocity of the earth, 0.25068 deg/min (taken as 

0.25 deg/min in this report) 

Subscripts : 

N nodal point 

n orbit number 

S any satellite position 

1 initial or burnout position 

2 selected position 

2e equivalent selected position in first orbit 

l-2e increment between initial or burnout position and equivalent 

position 

SATELLITE POSITION RELATIONSHIPS 


If drag or accelerating forces acting parallel to the path of the 
satellite are neglected, the initial conditions of speed Vp, radial 
distance rp, and elevation angle 7^ at final-stage burnout determine 
the orbit characteristics of a satellite. It is the inclination angle 
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of the orbital plane as determined by the azimuth angle ^ at the burn- 
out latitude however, which determines the range of latitude that 

covered by the satellite. It should, therefore, be possible for 
a given value of V]_, r^, and t ° determine an initial azimuth 

angle *^2. ' w ’^ich will place a satellite vehicle over a selected position 

on the rotating earth surface after a selected small number of orbits. 

A method for determining this initial azimuth angle Is presented in this 
report . 

The orbit relations obtained from elementary celestial mechanics 
for a body acted upon by a radial force of attraction proportional to 
the inverse square of the distance are first reviewed. Reference con- 
ditions for relating the position of the satellite In the orbital plane 
with the projected position on the earth surface are established. These 
results are then modified so that the latitude and longitude of the 
satellite trace on a rotating earth surface may be obtained. An expres- 
sion is derived which will establish the azimuth angle at burnout required 
to place a satellite over a selected earth position. Expressions for the 
incremental changes from the desired latitude and longitude due to errors 
in the azimuth angle and to the earth* s oblateness are presented in 
appendix A. 


Elliptical Orbit Relations 

From elementary orbital mechanics (ref. l), the characteristics of 
an orbit in terms of the initial conditions at burnout r-^, V^, and 7^ 

can be obtained from the following expressions: 

The equations for distance of satellite from earth center (see fig. l) 
are 


r 


P 

1 + e cos 0 


(la) 


a(l - e^) 
r = — - — 

1 + e cos 9 


(ib) 




r = a(l - e cos E) 


(lc) 
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where 


p r A 2 

— = — cos^7 

r l Vc 


Satellite velocity is given by 


v* - - - 

\r a 


Satellite velocity for a circular orbit at radius r-^ is 


2 *o r 


Angle in orbit plane between perigee and burnout point is given by 


tan = 


- 1 

n 


Eccentricity is 


e - - l' 


COS 9^1 


Satellite period in seconds is 


T = 2it 


R .//a 


Sn R 
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Time from perigee is expressed as 


t(e) = — (E - e sin E) (8) 

2rt 

Eccentric anomaly (fig. 1(b)) is given by 


E = 2 tan -1 




(9) 


Relations Between Orbit and Earth Systems of Coordinates 

In the present paper spherical-triangle relations are used for 
relating the satellite position in the orbital plane to its trace on the 
earth surface. By using the Inclination angle i of the orbit plane 
(which is numerically equal to the maximum latitude reached) and a ref- 
erence or nodal point N, it is possible to obtain the desired position 
relations between the equatorial and orbital planes. The angular motion 
of a satellite in an elliptic orbit is usually related to a perigee pas- 
sage j therefore, in the present paper time is referenced to the perigee 
point. The orientation of a satellite orbit and the projection of the 
orbit plane on a spherical earth are illustrated in figure 2. 

Latitude position .- The inclination angle i, the angle between the 
orbit plane and the equatorial plane (fig. 2(a)), is fixed at the instant 
of burnout by the latitude 0 1 and azimuth tp- If the equations for a 

right spherical triangle are applied to the triangle bounded by the equa- 
tor, the projected orbit plane, and the meridian or longitude through the 
position 0p, the pertinent relation for the situation pictured in fig- 
ure 2(b) is 


cos i = cos 0p sin \|rp (0° < \|rp < l 80 °) (10a) 


For a retrograde or westward launching, which occurs when the azimuth 
angle at burnout exceeds 180°, the relation is 


cos i 


-cos 0-^ sin ijr^ 


(l 80 ° < ijfp < 360 °) (10b) 


S 
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The position of the satellite in the plane of the orbit at the 
instant of burnout is specified by the angle 0p (eq. (5)) measured 
relative to perigee. The perigee point in turn is related to one of 
the two nodes at the intersection of the projected orbital plane and 
the equatorial plane by the single co. The relationship for the situa- 
tion shown in figure 2(b) is 


sin 0-, 

— — — = sin(co + 0p) (lla) 

sin i 

or 

tan 0-j 

-j r = tan(co + 0p) (lib) 

cos 

For determination of the argument of perigee <d, with 0]_, 0-^, and 

or i given, the inverse angle relations implied by equations (lla) 

and (lib) are ordinarily most conveniently resolved in favor of placing 
the angle (ju + 0^ in the first or fourth quadrants - that is, 

90 ° > (to + 0p) > - 90 °) thus, for reference purposes the node which is 

closest to the point of burnout is selected. 

Subsequent to burnout, the latitude of the satellite at any point 
in the orbit 0g is given by relations similar to equation (lla), that 
is 

sin 0g = sin i sin(co + 0s) (12) 

For use of equation (12) it is convenient to regard the* angle i as 
positive if the reference node is an ascending node, one where the equa- 
tor is crossed from south to north. Conversely, if the reference node 
is a descending node, i should be regarded as a negative angle. The 
combinations of launch direction and burnout latitude, which lead to 
positive or negative values of i, together with the range of the 
azimuth are summarized in table I. 

Longitude position .- The longitude of a point on the earth's surface 
is measured by the arc of the equator intercepted between the prime merid- 
ian at Greenwich and the meridian through the point. A convenient refer- 
ence for determining the longitude of a satellite is the node, the point 
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of intersection of the equatorial and orbital planes, at the instant of 
burnout. At a subsequent time, the nodal point will have shifted west- 
ward by the product of the earth's angular velocity eng and the orbiting 

time. 


As given in equations (8) and (9), time is regarded as a function 
of the orbital angle 0 and is conveniently computed from perigee, 
0=0°. In the present paper, orbiting time is measured from burnout, 

0 = 0]_. Time, therefore, at any instant in the first orbit may be given 
by t(0s) - t(0q) or, in general, by t(0s) - t(0i) + nT where T is 
the period of the satellite in minutes, and n is the number of com- 
pleted orbits. (An orbit is completed when the latitude at burnout 0^ 

is reached again in the same direction as at launch.) 

From figure 2(b), the longitude of the reference node (^N) re f 811(1 
the longitude of burnout are related by 

M r ef = - M N1 


where 


tan = sin 0^ tan llr^ 


( 1*0 


From figure 2(c), at time t(0s) 


tan AAjfg = cos i tan(o) + 0g) 


(15) 


and the longitude of the node is given by 


= ( a n) 


ref 




t(0s) - t(0i) + nT 


(16) 


Hence, the longitude of the satellite at any time is 

... _ . _ . . T 

A s = (^n) re f - (%&T + t(0 S ) - t(0i)] + AX NS 


(17) 
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EXPRESSION FOR PLACING A SATELLITE OVER A SELECTED 

POSITION 


The requirement that a satellite with burnout position 0q,Aq pass 
over a selected position 0 2 , ^2 after the completion of n orbits is 

equivalent to the requirement that, during the first orbit, the satellite 
pass over an equivalent position with latitude 0 2 the same as that of 

the selected position but with longitude A 2e displaced eastward from 
A 2 by an amount sufficient to compensate for the rotation of the earth 
during the n complete orbits, that is, by the polar hour angle najgT. 

(See fig. 3 (a)*) The longitude of this equivalent position is thus given 
by the relation 


A 2e = A 2 + nojgT 


(18) 


The azimuth angle which will place the satellite over the position 

0 2 ,A 2e at time t(0 2e ) will, therefore, n orbits later place it over 
the selected position 02>^2- 


The trace of the satellite position from the point 0q,Aq to the 
point 02>^2e is n °t a great circle, since during the time interval 
1(026) - t( 0q) when the satellite passes through the orbital angle 
e 2e “ 9 1> the earth rotates eastward through the polar hour angle 


cog t(02e) - t(0q)J . 


The difference in polar hour angle (which is a dif- 


ference in longitude) between the meridian Aq at t(0q) and the merid- 
ian through the point 0 2 ,A 2e at t(02 e ) will be 


^l-2e = x 2e “ ^1 + ^[^(^e) - t(0q)J (19) 

Since the latitudes 0q and 0 2 are known, the orbit angle 0 2e - 0q 
(fig. 3(h)), may be obtained from equation (19) and the relation 

cos(02e - 01 ) = sin 02 sin 0q + cos 02 cos 0q cos AAq_ 2e (20) 
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(in the application of equation (20) it may be noted that since the angles 
®2e “ ©l and A^l-2e &re both portions of intersecting great circles, 

©2e - 9p will be less than (or greater than) l80° when AAp_2 e is less 
than (or greater than) l80°.) 

In the use of equations (19) and (20) an iterative procedure is 
required, since the time t(02e) from perigee to the equivalent posi- 
tion is not known initially. A satisfactory first approximation is to 
assume that 


jt(e 2 e) - t(0ij] « ^L(A 2e - Ai) 


( 21 ) 


This value ; substituted into equation (19)* yields a first approximation 
for AA^_2e* This approximation substituted into equation (20) yields 

a first approximation for 02 e from which, by equation (8), a new value 
of t(02e) i s obtained. The iteration of these three steps is continued 
until the value of 02 e converges. For convenience these equations are 
listed in appendix B. 

The effects of the earth f s oblateness, discussed in appendix A, may 
be included as small modifications of the values of 02 anc ^ ^2e 

obtained in the first iteration. The equations are given in appendix B, 
and examples are shown in the following section. 

When convergence of the value of 02e has been attained, the 
az im uth angle at burnout ^ may now be calculated from the sine law 
(fig. 3(b)) by 


sin AA^_p e cos 02 
sin = 

sin(02e - ©l) 


and the inclination angle by equations (10): 


( 22 ) 


cos i = ±sin i|r^ cos 0-^ 


vA C\) 00 OA 


11 




The argument of perigee a o may he calculated from the following 
equation which is equation (lla): 


sin(o> + 0q) 


sin 0-^ 
sin i 


The position equations given in the previous section and listed in 
appendix B may he used to check final values and to compute the complete 
orhit path. 

In orhit calculations, such as those described in this paper, a 
number of occasions arise when ambiguities of the antitrigonometric 
functions must he resolved in order to place angles in the proper quad- 
rants. The following conventions have been used in this paper: 

Longitude A: range, 0° to 360°; positive, counterclockwise, east- 

ward; negative, clockwise, westward 

Azimuth >|r^: range, 0° to 3^0°; clockwise from north; orbits exhibit 

normal rotation when 0° < \(r^ < l80°; orbits exhibit retrograde motion 
when 180° < \|r^ < 3^0° 

Orbit angle 0: range, 0° to 360°; positive for normal rotation; 

negative for retrograde motion 

Inclination angle i: range, 0° to 90°; measured from equator; 

positive at ascending node; negative at descending node 

Nodal point N: sometimes specified as ascending node; reference 

node utilized herein is selected on basis of proximity to burnout point 


EXAMPLE ORBIT CALCULATIONS 


In order to Illustrate the application of the developed expressions, 
a temporary satellite orbit has been selected just outside the earth’s 
atmosphere . Examples are given of a satellite launched eastward and one 
launched westward, each passing over a selected position after having 
completed three orbits. The two directions of launch Illustrate that the 
procedure is not limited by the direction of launch. The positions were 
considered favorable positions from which to initiate a reentry leading 
to a selected recovery area. Incremental changes from the desired lati- 
tude and longitude due to the earth's oblateness, as discussed in 
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appendix A, are included in the procedure for obtaining the azimuth angles. 

The results for both cases are then compared with the results obtained 
by the method of reference 2 of a three-degree-of- freedom calculation 
for an oblate rotating earth performed on an IBM 70^ electronic data 
processing machine. The data from this calculation are referred to here- 
inafter as three-degree-of-freedom data. Changes from the selected lati- 
tude and longitude, resulting from incremental changes from the burnout 
condition \|r-^ and 0p are also analyzed. 

The orbit characteristics of a satellite vehicle launched into an 
orbit 120 nautical miles above the earth's surface with an initial veloc- 
ity V]_ = 1.010V C and an elevation angle y = 0.50° are presented in 

table II. The satellite radial distance r and time t calculated by 
equations (l), (8), and (9) are plotted as a function of 8 in figure 4. 9 

With the three parameters - speed, radial distance, and elevation angle - 
given, only the azimuth angle at burnout \|r-]_ remains to be calculated 

to establish the orbit completely. 

Eastward Launch, Case A 

The satellite is to be launched eastward from a burnout position 
01=28.50° N. ,A^=279*^5° E. and is to pass over a position 
02=3^.00° N. ,A2=24l.00° E. after completing three orbits. The steps of 

the procedure for determining the azimuth angle at burnout are outlined 
and illustrated in appendix B. 

The azimuth angle \tr-j_ is determined on the basis of an equivalent 
position 02 >^2e which must be passed in the first orbit. After the 

first iteration. It is desired to incorporate the oblateness effects as 
corrections to the desired position and to 02>^2e* For the orbiting 

time 3T + t(02e) - t(0i) « 281.615 minutes, there results a change in 
latitude of 0.086° and a change in longitude of 1.028°. These values 
are subtracted from the original latitude and longitude to give a cor- 
rected position 02=33 *91^° N. ,A2=239-972° E. and hence a corrected 

A 2 e=308.66l° E. With these substitutions into equations ( 19 ) and (20), 
the iteration for ©2e is continued and final values of the desired 
quantities are obtained as shown in appendix B. 

By using the given initial conditions and +1 calculated by the 

method of this paper, an orbit was calculated by a three-degree-of- 
freedom method for an oblate rotating earth (ref. 2). A trace of the 
resulting orbit is presented in figure 5(a). The method of the present 
paper indicated that the selected position would be reached in 


(-^ 0J00 
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28l.5l6 minutes. The three-degree-of-freedom data indicated that at 
t = 281.516 minutes the satellite is located at 0=55-98° N. ,A=2l0.62° E., 
which is l 8.6 nautical miles west of the desired position, and would reach 
the desired position at t = 281.598 minutes. This agreement is con- 
sidered good. 

With the calculated azimuth and inclination angles, errors in the 
desired latitude and longitude due to insertion errors in and 0 -^ 

are examined. The results for case A are presented in figure 6 (a). The 
variations of Ai, A 02 , and AA 2 with A\|r^ and A0]_ are considered 

to he small, the slopes being dependent on equations (Al) and (A5). 


Westward Launch, Case B 

The same iteration procedure is used for the westward launch as 
used for the eastward launch with pertinent sign changes . As shown in 
appendix B, first approximations are made and resulting values are used 
to compute oblateness effects. These corrections are applied to 02 

and A 2 , and final values are then derived. 


As in case A, an orbit was calculated by a three-degree-of-freedom 
method by using the given initial conditions and ^ calculated by the 

method shown in appendix B. 

The trace of the resulting orbit is shown in figure 5(b). The method 
of this paper indicated that the desired position would be reached in 
555-608 minutes. The three-degree-of-freedom data indicated that at 
t = 355.608 minutes the satellite is located at 0=20.63° N.,A=320.67° E., 
a point 52.2 nautical miles west of the desired position. 

The variations of Ai, A0 2 , and AA2 with Af x and A0-, for 

case B as derived in equations (Al) to (A3) are presented in figure 6 (b). 
The results for both cases indicate that the error in the desired lati- 
tude and longitude is a function of the orbital angle cd + 02 , the azimuth 
angle at burnout and the burnout latitude 0 1 . 


CONCLUDING REMARKS 


It is possible by the use of a simple iteration procedure to estab- 
lish the azimuth angle at burnout on the basis of a selected latitude and 
longitude and a prescribed number of complete orbits for a given satellite 
period and orbit eccentricity for either an eastward or westward launch. 



Incremental changes from the desired latitude and longitude positions due 
to the perigee motion and rotation of the orbital plane can conveniently 
be t ak en into account in the procedure. The magnitude of the changes for 
a given orbit is a function of the time in orbit and the inclination 
angle of the orbital plane. A comparison of the results for both the 
eastward and westward launches with the results obtained from a three- 
degree -of- freedom calculation for an oblate rotating earth shows very 
good agreement. This indicates that there are no apparent limitations 
as to launch direction or selected passage point. The error in the 
desired latitude and longitude is a function of the orbital angle, the 
azimuth angle at burnout, and the burnout latitude. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., November 27, 1959- 
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APPENDIX A 

V 

INCREMENTAL CHANGES FROM DESIRED POSITIONS 


Incremental changes In the azimuth angle ^ and the resulting 
Influence on the desired satellite passage point are P resent ^d 

along with incremental changes due to the effect of the earth’s oblate- 
ness. Changes In the desired burnout conditions r-^ and 7^ 

Introduce changes in the orbit elements which are beyond the scope of 
this report and will therefore not be considered. 

The error in the inclination angle Ai obtained by differentiating 
equation (lOa) with respect to ^ and is 

cos 0L cos iK sin >| a sin $ ^ 

A1 sin 1 A *1 + ^n— 1 4*1 < A1 > 


The error in the desired satellite latitude obtained by differentiating 
equation (10a), making use of equation (Al), and neglecting A0 2e is 


, sin(o) + 0 2p ) 

^ " 1 cos T~ tan 1 ( Sln *1 Sln *1 - COS *1 C0S *1 A *l) 


+ 


sin i cos(o> + 
cos 


Ad 


(A2) 


Similarly, the error in the desired longitude obtained by substituting 
equation (15) into equation (17) and differentiating, making use of 
equation (Al), and neglecting A0g e is 


aa 2 = 


tan(cu + 0 p ) 

-1 — ^ e/ (cos q , , cos t-, Avh - sin \|r-, sin 0, A0-, ) 

1 + cos^i tan 2 ^o) + 0 2g ^' X ^ 


cos i sec 2 ^oo + 0 2e ^ 

1 + cos 2 i tan 2 (u> + 0 2e ) 


■ An 


(A3) 
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The errors in (j > 2 and A 2 as presented by equations (A2) and (A3) are 
thus made up of errors in the insertion conditions aa d and 

errors due to perigee motion Ad. These errors are as follows: 

*2 - (^a) le + (*2)^ 

and 

AA 2 = (AA 2 ) ie + (^2) pm 


where the subscripts ie and pm represent insertion errors and 
perigee motion (oblateness effects), respectively, and 



sin i cos (ad + 0 O } 

= ± \ An 

cos 0 O 


(A4) 


and 



cos i sec^(a> + 0 2e ) 

OO V 

1 + cos^i tan ^o) + 0 2e j 


(A5) 


The rotation of the perigee point Ad is defined as 


An = — At 
dt 


(A6) 


where the approximate mean rate of the argument of perigee motion or 
rotation of the major axis in deg/min (ref. 3) is 


do> 

dt 


3/2. 


* 3.1*722 x 10-5(1) (I) (5 


o 

cos i - 


1 ) 


(A7) 


The rotation of the major axis is in the same direction as the satellite 
if i < 63.4 0 , in the opposite direction if i > 63-4° and is zero when 
i = 63.4 0 . The other main effect of the earth's oblateness is the rota- 
tion of the orbital plane A ft defined as 

Aft = — At (A8) 
dt 


h^l OJ co o\ 
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The orbital plane rotates about the earth's axis in the opposite direc- 
tion to the earth's rotation and the approximate mean rate of rotation 
of the orbital plane in deg/min (ref. 3) is 


ff = - 6 . 9 W+ X 10 _5 (!) COS i (A9) 


The rotation of the orbital plane affects only the incremental longitude 
so that a Aft term must be included in equation (A3) to give the total 
oblateness effect on the longitude of the desired passage point. The 
magnitudes of Ajo and Aft are obtained by multiplying equations (A7) 
and (A9) by the orbiting time nT + t ( 0 ) g - t(0)-^ expressed in minutes. 

Changes from the desired satellite passage point due to oblateness 
effect can be taken into account in the iteration procedure previously 
discussed. The desired passage point becomes 


02 ~ (A0p) y - (AAp) 

* V ^ 'total d v total 


These values are substituted into equation (19) before a second itera- 
tion is performed. 
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appendix b 

EQUATIONS FOR CALCULATING AZIMUTH ANGLE, OBLATENESS EFFECTS, 

AND SATE LLI TE POSITION AND ILLUSTRATIVE EXAMPLES 

Azimuth Angle 

The calculation of the azimuth angle at burnout involves the use 
of equations (21), (l8), (19), (20), (8), (22), and (10) as follows: 

[t( 0 2e) " t(8 1 )] * ^( A 2e " A l) 

M l-2e = ( A 2 + KD E T ) + “E [ t ( 0 2e) ’^l)] “ A 1 


COS 


- 0^ = sin $ 2 s ^ n $1 + cos ^2 cos ^1 cos ^l-2e 


t( s 2e) ' S (E ' 6 Sln E> 

sin A 2e cos $ 2 

310 *1 = ^(4 - »J~ 

cos i = ±sin \|/^ cos ^ 


Oblateness Effects 

The oblateness effects are calculated through the use of the fol- 
lowing equations: 

Ajd = 3.4722 X 10 " 5 (§) (|) cos 2 i - l)jnT + t(0 2e ) - ^(ejJ 

sin i cos^cu + 8ge) 


K4 = 1 


COS tyr 




►-» oj 00 cr\ 


rvi m cr\ 
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Aii = 


2 3/2 

-6.9410. x 10' 3 (|) (|) cos i[nT + t(0 2e ) - t(6 1 )] 


(AA 2 ) 


pm 


cos i sec^a) +©2e) 

2 " p 

1 + cos i tan^co + 0g e ) 


An + A Q 


These equations make use of equations (A5) to (A9) . 


Satellite Position 

The calculation of the satellite position involves the use of 
equations (12), (17), (lla), (l4), (13), and (15) . For convenience, 
these equations are repeated as follows: 


sin = sin i sin(o> + 0g) 
A S = ( X N) ref + AA ns - ^[nT + t(9g) 


sin^co + 0-jJ = 


sm 


h 


s In i 



tan = sin (jf^ tan \|r^ 


( X N) ref " A 1 - M N1 


tan AAjjg = cos I tanfa) + 


(w + 0g) 


Illustrative Examples of Orbit Calculations 

In the following sections examples of the orbit calculations for 
cases A (eastward launch) and B (westward launch) are given to demon- 
strate the application of the method. 
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Case A, eastward launch .- For case A, eastward launch, the given 
launch position is 02 = 28 . 50° N. ,>^=279 -^5° E.j the selected position 

is 0 2 =54-.OO° N. ,A 2 =24l.00° E.j the number of orbital passes n is 3j 
and the equivalent selected longitude A 2g is 309 .689° E. The values 
of the first approximations are 

[t( e 2e) - t ( e l)] “ I' 6 ” 

AA-L^g = 32.162 

0 2e = 51*89 
t( 0 2e) = 12.702 
t = 70.468 
i = 34.081 

The values found in the calculation of the oblateness effects are 

Ad = 1.965 
A0 2 = 0.086 
= -1.340 
AA 2 = 1.028 

Corrected 0 2 = 33*914° N. Corrected A 2 = 239*972° E. 

The iteration procedure, as explained in the text of this paper, yields 
the following values: 
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Parameter 

2d iteration 

3d iteration 

4th iteration 

t ( 0 2e) “ t ( 9 l) 

6.860 

6.603 

6.594 

^l- 2 e 

30 . 926 

30.862 

30.859 

9 2e 

50.866 

50.809 

50.806 

1 *(®2e) 

12.445 

12.436 

12.436 


The final values are as follows: 


^ = 70-541 
i = 34.043 

(a k ) = 225.971 

” ref 

03 = 34.497 

Case B, westward launch .- For case B, westward launch, the given 
launch position is 0 ] _=34.OOO N. ,7^=241.00° E.; the selected position 

is ^2 = 20.00° N.,A 2 =320.00 o E.j the number of orbital passes n is 3j 
and the equivalent selected longitude A 2e is 28.689° E. The values 
of the first approximations are 

t( 0 2e “ 0 i) = 54.012 

AA l-2e = -!98.808 

9 2e = “ 260. 863 

t( 0 2e) = 66.996 

^ = 201.203 

i = -72.555 


The values found in the calculation of the oblateness effects are 
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Ad = -0.531 

A$ 2 = 0.505 

AQ = -0.578 
AA 2 = -0.759 


Corrected <f > 2 = 19-497° N. Corrected X 2 = 320.759° E. 
The iteration procedure yields the following values: 


Parameter 

2d iteration 

3d iteration 

4th iteration 

t( 0 2e) “ t ( Q l) 

61.154 

60.848 

60.853 

^1-26 

-196.263 

-196.340 

-196.339 

0 2e 

-259.664 

-259.685 

-259-684 


66.690 

66.695 

66.695 


The final values are as follows: 

if 1 = 198.721 
i = -74.569 
(*») ref - 2 5°*269 


a) » 59.427 


^ CVI CO CT\ 



OJ oo o\ 


23 


REFERENCES 


1. Moulton, Forest Ray: An Introduction to Celestial Mechanics. Second 

rev. ed.. The Macmillan Co., c.1914. 

2. Nielsen, Jack N. , Goodwin, Frederick K., and Mersman, William A.: 

Three-Dimensional Orbits of Earth Satellites, Including Effects of 
Earth Oblateness and Atmospheric Rotation. NASA MEMO 12-4-58A, 
1958. 

3. King-Hele, D. G., and Gilmore, D. M. C.: The Effect of the Earth's 

Oblateness on the Orbit of a Near Satellite. Tech. Note 
No. G.W.475, British R.A.E., Oct. 1957- 



TABLE I 


MAGNITUDE OF AZIMUTH ANGLE AND SIGN OF INCLINATION ANGLE i 

FOR A GIVEN BURNOUT LATITUDE AND LAUNCH DIRECTION 
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TABLE II 

ORBIT CHARACTERISTICS USED FOR SAMPLE CALCULATIONS 

Satellite velocity at burnout, V-^, ft/sec 25761-345 

Distance of satellite from center of earth at 

burnout, r^, ft 21, 637 >933 

Elevation angle, deg 0-5 

Velocity for circular orbit, V c , ft/sec 25506.28 

p/^ 1.020022269 

Semimajor axis, a, ft ..... 22081775*58 

Angle between perigee and burnout, 0-^, deg ±23.969 

Eccentricity, e 0.0219118 

Period, T, min 91*585 

Time from perigee to burnout, min 5*842 
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(b) Relationship between orbit angles and earth positions at 

time t(02e)* 


Figure J>.~ Relation between selected position and equivalent position. 
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upn f (e)* < 883Ta©d od.zj 


Figure 4.- Satellite time and distance as a function of 0. r apogee ~ 22565628 ft 

r perigee = 21597924 ft; T = 91-585 min. 
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(b) Case B, westward launch* 


Figure 5*“ Concluded. 







